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ABSTRACT 

This paper presents 13 color CCD intermediate-band spectrophotometry of 
a field centered on the open cluster M48 (NGC 2548), from 400nm to nearly 
lOOOnm, taken with Beijing-Arizona- Taiwan-Connecticut (BATC) Multi-Color 
Survey Photometric system. The fundamental parameters of this cluster are 
derived with a new method which based on the comparison between the spectral 
energy distributions (SEDs) of cluster stars and the theoretical SEDs of Padova 
models. We find that the best-fitting age of M48 is 0.32 Gyr, the distance is 780 
pc and a reddening E(B — V) = 0.04 with a solar metallicity Z = 0.019. 

Subject headings: Galaxy: open clusters and associations: individual: M48 (NGC 
2548) — Stars: general — Hertzsprung-Russell diagrams 
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1. Introduction 

Open clusters have long been recognized as important tools in the study of the Galactic 
disk. Their value lies in the improvement in accuracy for distance determination, metal 
content, reddening and age produced by the collective stellar sample which shares these 
properties. So open clusters are excellent tracers of the abundance gradient along the 
Galactic disk as well as many other important disk properties, such as the age-met allicity 
relation, abundance gradient evolution, disk age, and so on (Friel 1995; Twarog, et al. 1997; 
Chen et al. 2003). 

The open cluster M48, also known as NGC 2548 is a quite conspicuous object and 
should be a naked-eye object under good weather conditions. It was firstly studied by 
Peshch (1961) who made UBV photoelectric photometry of 37 stars, from which the 
cluster turned out to have a reddening E(B — V) — 0.04 ± 0.05 and a distance 630 pc. 
DDO photoelectric photometry of 5 red giants was obtained by Claria (1985). From 4 
members of these red giants, he conclude that this cluster has a metallicity [Fe/H]=0.1, 
E(B -V) = 0.06 ± 0.02 and a distance of 530 pc. Mermilliod (1981) derived an age of 0.30 
Gyr based on a synthetic composite color-magnitude diagram. Wu et al. (2002) determined 
absolute proper motions and membership probabilities for 501 stars in the field of M48. 
More recently, Rider et al. (2004) obtained a new result of M48 taken in the u'g'r'i'z SDSS 
filter system. They find that a distance of 700 pc, an age of 0.40 Gyr and a metallicity of 
[Fe/H]=0.0 can fit their data best. 

In this paper we present a new photometric result of M48 taken with BATC Multi-Color 
Survey Photometric system. The BATC filter system consists of 15 filters of band-widths 
150 - 350A that cover the wavelength range 3300A - 10000A, which avoids strong and 
variable sky emission lines (Fan et al. 1996). As the first object in BATC survey, the 
old open cluster M67 has been studied based on color-magnitude diagram (CMD)(Fan et 
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al. 1996). Using the BATC filter system, Chen et al. (2000) studied the globular cluster 
NGC 288 by comparing SEDs of bright stars with Kurucz models. The estimated effective 
temperatures and average value of [Fe/H] for these stars are consistent with spectroscopic 
determinations. Based on the BATC survey observations, the main aim of this study is 
to determine simultaneously the fundamental parameters of M48, such as age, distance, 
metallicity and reddening, by comparing observational SEDs of cluster stars with theoretical 
stellar evolutionary models. 

The observations and reduction of the M48 data are described in Sec. 2. In Sec. 3, we 
derive fundamental parameters of M48. Conclusions and summary are presented in Sec. 4. 

2. Observation and Data Reduction 
2.1. Observation 

The observations are done with BATC photometric system at Xinglong Station of 
National Astronomical Observatories, Chinese Academy of Sciences (NAOC). The 60/90 
cm f/3 Schmidt telescope is used with a Ford Aerospace 2048 x 2048 CCD camera at its 
main focus. The field of view of the CCD is 58' x 58' with a plate scale of 1.7 arcsec per 
pixel. 

The filter system of BATC project is defined by 15 intermediate-band filters which 
are designed specifically to avoid most of the known bright and variable night sky emission 
lines. The definition of magnitude for the BATC survey is in the AB V system which is a 
monochromatic F u system first introduced by Oke & Gunn (1983): 

m batc = ~ 2 - 5 lo S F » ~ 48 - 60 i 1 ) 
where F v is the appropriately averaged monochromatic flux (measured in unit of 
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erg s 1 cm 2 Hz 1 ) at the effective wavelength of the specific pass-band (Fukugita et al. 
1996). In BATC system, the F v is defined as (Yan et al. 2000): 

_ fd(\ogv)f v R v 
* v !&{\ogv)R v [Z) 

which ties directly the magnitude to input flux. The system response R u actually used to 
relate the spectrum energy distribution of the source f u and F u , includes only the filter 
transmissions. Other effects, such as the quantum efficiency of the CCD, the response of the 
telescope's optics, the transmission of atmosphere, etc., are ignored. This makes the BATC 
system filter-defined, and the effective wavelengths are affected only at the ±6A level after 
taking CCD quantum efficiency and aluminum reflection into account (Yan et al. 2000). 

The flux calibration of the BATC photometric system is defined by four spectropho- 
tometric standard stars of Oke k Gunn (1983): HD 19445, HD 84937, BD+26°2606 and 
BD+ 17°4708. The fluxes of the 4 standard stars have been re-calibrated by Fukugita et al. 
(1996). Their magnitudes in BATC system have been slightly corrected recently by cross 
checking with the data obtained on a number of photometric nights (Zhou et al. 2001). 

In the nights judged photometric by the observers, the standard stars were observed 
between air-masses 1.0 and 2.0 for each filter band. The observing procedure of survey 
program field and photometry are described in detail in Zhou et al. (2001) and Yan et al. 
(2000). 

Because of the very low quantum efficiency of the thick CCD used in the bluest filters, 
two BATC filters a and b are not used in the observation of M48 field. In Table 1, for each 
BATC filter, we list the corresponding effective wavelength, FWHM and exposure time. 
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2.2. Data Reduction 

Preliminary reductions of the CCD frames, including bias subtraction and field 
flattening, were carried out with an automatic data reduction procedure called Pipline I 
which has been developed as a standard for the BATC survey in NAOC (Fan et al. 1996). 
The astrometric plate solution is obtained by a priori knowing the approximate plate center 
position, and then using this information to register the brighter stars in each frame with 
the Guide Star Catalog (GSC) coordinate system (Jenkner et al. 1990). 

A Pipeline II program which based on the DAOPHOT II stellar photometric reduction 
package of Stetson (Stetson 1987) was used to measure the instrumental magnitudes of 
point sources in BATC CCD frames. The Pipeline II reduction procedure was performed 
on each single CCD frame to get the PSF magnitude of each point source. The instrumental 
magnitudes were then calibrated to BATC standard system (Zhou et al. 2003). The average 
calibration error of each filter is less than 0.02 mag. The other sources of photometric error 
include photon statistics of star and sky, readout noise, random and systematic errors from 
bias subtraction and flat-fielding, and the PSF fitting, are considered in Pipeline II and the 
total estimated errors of each star are given in the final catalog (Zhou et al. 2003). Stars 
which are detected in at least 3 filters are included in the final catalog. 

3. Fundamental Parameters Derived from SEDs 

CMD is the main tool to derive the fundamental parameters of star cluster. In the 
BATC photometric system, 15 filters can be used to form various CMDs. However, it is 
difficult to derive a consistent result from various different CMDs. In the other hand, 
the large number of bands available in our data set provides a sort of low resolution 
spectroscopy, which defines the SED of each star quite well. So it is possible to derive the 
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fundamental parameters of star cluster by fitting the SEDs of member stars with theoretical 
stellar evolutionary models. 

3.1. Fitting Procedure and Results 

The observed SED of a star is determined by intrinsic (mass, age, metallicity) and 
extrinsic (distance and reddening) parameters. The member stars in a star clusters share the 
same age, metallicity, reddening and distance. The idea of our new method is to compare 
the observed SEDs of member stars with theoretical models to obtain a combination of 
best-fitting cluster parameters. 

Our fitting procedure can be separated into two steps. First, we calculate the deviation 
between observed and theoretical SEDs of each member star with different sets of cluster 
parameters including age, metallicity, distance and reddening: 

m x [E(B -V),r] = m Xoba + 5 - 5 logr - R x x E(B - V) (3) 

where m Xobs is the observed magnitude of a cluster member in filter x, m x [E(B — V),r] is 
the absolute magnitude corrected by distance r and reddening E(B — V), R x is extinction 
coefficient which transform the E(B — V) to the BATC filter system and derived by Chen 
(2000) based on the procedure in Appendix B of Schlegel et al. (1998). 

The criteria of fit for a cluster member s between observed SED and theoretical model 
is defined by: 

(.M.^-viH ^B E (4) 

where M x is the theoretical magnitude of a star in filter x, computed from an chosen 
theoretical isochrone models with age log(t), metallicity Z and mass m, n is the number of 
filters used in the fit. 



- 8- 



For any set of cluster parameters combination, we can find the smallest value of ( s viz. 
MIN^ s [\og(t),z,r,E(B-v)] f° r the member star s and give the best-matched theoretical mass of 
that star. Then we can get the total deviation of all member stars of cluster under this set 
of fundamental parameters. 

C [log( t) . Z, r, E(B - V)l = E MIN <- (5) 

where N is the number of cluster members used in the fit. For a set of fundamental 
parameters that best match the observed and theoretical SEDs of cluster members will 
yield the minimum of £• 

As our request, Dr. L. Girardi has kindly calculated isochrones of our filter system 
using the known BATC filter transmission curves and their Padova stellar evolutionary 
models (Girardi et al. 2000, 2002). The Padova isochrone sets are computed with updated 
opacities and equations of state, and a moderate amount of convective overshoot. 

The results of Wu et al. (2002) proper motion and membership study of M48 are used 
to determine members in this cluster. Stars with membership probabilities greater than 0.7 
are considered to be members (Wu et al. 2002). All stars considered as members based on 
their proper motions were used in our fitting. 

In our fitting procedure, the distances was chosen from 600 pc to 900 pc in intervals of 
10 pc, E(B — V) from 0.00 to 0.10 in intervals of 0.01. Theoretical isochrone models with 
metallicity Z = 0.08, 0.019, 0.030 and age log(t) from 8.0 to 9.0 in intervals of 0.05 were 
chosen. 

We find that, with a distance of 780 pc, reddening E(B — V) — 0.04, the theoretical 
model with age of 0.32 Gyr and metallicity Z = 0.019 give the smallest value of ( and best 
fit the observed SEDs. In Figure 1, we plot the best-fitting results for some member stars. 
The mass of each object in Padova models is labelled on the right of each corresponding 
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curve in the unit of solar mass M Q . In the top panel of Figure 1, SEDs of 11 main sequence 
(MS) stars with mass from 1.5662 M Q to 3.1878 M & are plotted. In the bottom panel of 
Figure 1, the SED of a red giant star is plotted. We can see that the observed SEDs of both 
MS stars and red giant stars with the derived best-fitting parameters, can fit the theoretical 
ones very well. 

In Figure 2, we plot four representative CMDs from our data: (c — p) vs c gives us the 
widest pass-band colors; (c — e) vs c gives us the cleanest CMD; and (/ — i) vs / give us the 
deepest CMD. All of the CMDs have a well defined MS and MS turnoff point. All stars in 
the field of M48 are plotted in each diagram. Stars with known membership probabilities 
greater than 0.7 determined from the proper motion study of Wu et al. (2002), are given 
distinct plotting symbols. In each CMD of Figure 2, parameters derived from SED-fitting 
are adopted. The theoretical isochrone with the parameters derived by SED-fitting can fit 
star's distribution in each CMD very well. 

The uncertainties in our derived best-fitting fundamental parameters can be determined 
from the observational photometric errors. For each member of cluster that used in previous 
fitting procedure, a new magnitude in each BATC filter was generated in a Monte Carlo 
fashion by adding Gaussian deviates to the observed magnitude. The standard deviations 
of the deviates in each BATC filter were taken to be 0.05 mag which is the maximum 
photometric error in most filters of used sample stars. All of these stars with new artificial 
magnitudes were then fitted by the same procedure as previous section. We repeated this 
simulation 100 times. In the end, we find that the uncertainty caused by photometric errors 
in distance is ±10 pc, in reddening E(B — V) is ±0.01, in age log(t) is ±0.05. The obtained 
metallicity keeps the same value. So the effect of photometric errors on the best-fitting 
results is very small. At the same time, there still remains the possibility that there could 
be some sort of systematic calibration problem between the photometry and the models, or 
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systematic errors in this particular set of photometry, although such errors are small, and 
not necessarily zero. 

In Table 2 we list the fundamental parameters of M48 derived from this and previous 
works. The derived age of this work is consistent with those derived by previous works, the 
largest difference is 0.08 Gyr (Rider et al. 2004). Our derived distance is very close to these 
results derived in resent years (Dias et al. 2002; Rider et al. 2004). Our derived metellicity 
is consistent with that of Rider et al. (2004) and close to the value derived by Claria (1985). 
The reddening E(B — V) derived in this work is same as that of Peshch (1961) and very 
close to those derived by previous works. 



4. Conclusions 

In this paper, we present and discuss new BATC multicolor photometry results for 
the intermediate age open cluster, M48. Comparing the observed SEDs of cluster member 
stars with the theoretical SEDs of Padova models, we find a set best-fitting fundamental 
parameters for this cluster: an age of 0.32 Gyr, a distance of 780 pc, a metallicity of 
Z = 0.019 with a reddening E(B - V) = 0.04. This SED-fitting result can also fit CMDs 
formed from our data very well. Our derived values are also consistent with those derived 
by previous authors (Mermilliod 1981; Rider et al. 2004). So, we can say that using the 
Padova theoretical isochrones, we can effectively fit our data obtained based on BATC filter 
system with theoretical isochrones and SEDs to extract useful cluster information. 
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suggestions that improved this paper. We wish to thank Dr. Leo Girardi for his assistance 
in the theoretical models. Z.Y.W acknowledge Dr. Wen-Ping Chen for many useful 
suggestions. This work has been supported by the Chinese National Key Basic Research 
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Fig. 1. — SEDs of M48 members of BATC data compared with the Padova models. The 
open squares are the observations plotted at the effective wavelengths of the BATC filters. 
The line is the best-fitting theoretical SED from Girardi et al. (2000). For each line, the 
corresponding mass is labelled on the right in the unit of solar mass M . 
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Fig. 2. — The color-magnitude diagrams for M48. Stars with membership probabilities 
greater than 0.7 considered as members (Wu et al. 2002) are plotted as squares. A distance 
of 780 pc, E(B — V) = 0.04 and a metallicity of Z = 0.019 are adopted, a Padova isochrone 
with age of 0.32 Gyr is over-plotted in each CMD. 
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Table 1. Parameters of the 13 BATC filters and statistics of observations 



No. 


Filter 


A eff (A) 


FWHM (A) 


Exposure time (min) 


1 


c 


4180.0 


310.0 


38 


2 


d 


4532.0 


330.0 


183 


3 


e 


4916.0 


370.0 


98 


4 


f 


5258.0 


340.0 


84 


5 


9 


5785.0 


290.0 


64 


6 


h 


6069.0 


310.0 


41 


7 


i 


6646.0 


490.0 


44 


8 


3 


7055.0 


240.0 


11 


9 


k 


7545.0 


190.0 


21 


10 


m 


8020.0 


260.0 


76 


11 


n 


8483.0 


170.0 


88 


12 


o 


9180.0 


250.0 


148 


13 


P 


9736.0 


280.0 


134 
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Table 2. Fundamental parameters of M48 



Author 


Distance 


Age 


Reddening 


Metallicity 




pc 


Gyr 


E(B - V) 


[Fe/H] 


This work 


780 


0.32 


0.04 


0.0 


Rider et al. (2004) 


700 


0.40 


0.03 


0.0 


Dias et al. (2002) 


770 


0.36 


0.03 


0.08 


Claria (1985) 


530 




0.06 


0.1 


Peshch (1961) 


630 




0.04 




Mermilliod (1981) 




0.30 







